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Abstract 

The content of Einstein's theory of gravitation is encoded in the prop- 
erties of the solutions to his field equations. There has been obtained a 
wealth of information about these solutions in the ninety years the theory 
has been around. It led to the prediction and the observation of physi- 
cal phenomena which confirm the important role of general relativity in 
physics. The understanding of the domain of highly dynamical, strong 
field configurations is, however, still quite limited. The gravitational wave 
experiments are likely to provide soon observational data on phenomena 
which are not accessible by other means. Further theoretical progress will 
require, however, new methods for the analysis and the numerical cal- 
culation of the solutions to Einstein's field equations on large scales and 
under general assumptions. We discuss some of the problems involved, 
describe the status of the field and recent results, and point out some 
open problems. 



1 Introduction 



The fascinating reports on the exciting new theories, which propose to unite 
our present and all future ideas about space-time, gravitation, and quantum 
physics into one coherent scheme from which general relativity will be derived 
in the end as a particular limit, may suggest that Einstein's classical general 
relativity is essentially understood. This impression is easily corroborated by 
the amazing successes of general relativity. Einstein concludes his synopsis from 
1916 of the theory of general relativity ([32]) with three predictions: the red 
shift, the bending of light rays, and the precession of the perihelion in planetary 
motion. In view of the present observational facts the revolutionary character 

•Extended version of a talk which was to be delivered at the DPG Friihjahrstagung in 
Berlin, 5 March 2005. 
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these implications may have had at the time falls into oblivion. The red shift can 
be measured to day in terrestial experiments f |lllp . the bending of hght rays 
is used under the heading 'gravitational lensing' as an effective astrophysical 
tool and various double star systems provide laboratories for relativistic 

gravity which let the minute advance of the periastron of Mercury of 42.98 
arcsec/century fade into insignificance: in the case of the binary pulsar PSR 
1913 + 16 the observed advance of periastron is 4.226595(5) deg/yr f |136j ') and 
for the recently discovered double pulsar J0737-3039 it even amounts to 16.9 
deg/yr {^). 

In the articles |S7|, Einstein discusses 'gravitational radiation' and de- 
rives the famous quadrupole formula. Doubts have been raised subsequently 
whether the notion of gravitational radiation referred to a real physical phe- 
nomenon (cf. |104| ). but again the prediction has been confirmed convincingly. 
Using Einstein's quadrupole formula to calculate the rate of period decrease of 
the system PSR 1913 + 16 due to its emission of gravitational radiation, one 
obtains, after taking into account certain small corrections, a curve which shows 
an uncanny agreement (to within about 0.2 percent) with the data gathered over 
the last 30 years C |136j ). 

The global studies of general relativity starting in the second half of the last 
century, led, together with unexpected observations, to concepts which went far 
beyond what had been envisioned by Einstein. Among those the notion of a 
black hole, a pure space-time structure which has no place in special relativity, 
is certainly the most remarkable one. Though the derivation of detailed obser- 
vational information still poses difficulties, the present situation suggests that 
black holes have to be accepted as part of our reality (|SZ|j jll4j l. More could 
be said in support of general relativity but I shall leave it at that. 

The overwhelming success of general relativity alluded to above may suggest 
a clear and simple answer to the question posed in the title of the present article. 
Also, we understand, of course, what it means that the geometry of the world 
is modelled by a Lorentz metric g^j,^ on a 4-dimensional manifold M and that 
this structure is governed by Einstein's field equations 

R^,^ ^ K{Tf,,, gf,^) + Xg^^, (1.1) 

with cosmological constant A, together with equations for matter fields which 
define the energy-momentum tensor T^i, and its trace T. In fact, the predictions 
referred to above have been derived from these requirements. There is, however, 
still a large and potentially most important part of the theory we do not have 
access to, neither mathematically, nor theoretically, nor observationally. 

On the observational side this situation may change soon. With the gravita- 
tional wave experiments presently becoming operational, we may well enter one 
of the most important eras of experimental relativity. The mathematical or the- 
oretical side seems to be lacking behind, however. In spite of a few general and 
a rich collection of specific results our qualitative and quantitative understand- 
ing of highly dynamical processes, strong field situations, and Einstein evolution 
over long time scales under general assumptions is still quite limited. 
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This may sound more pessimistic than it should. The field has seen a sub- 
stantial progress in recent years and for young researchers there is a unique 
opportunity to contribute to the theoretical investigation of an unexplored do- 
main of fundamental physics. Confronted with the gravitational wave measure- 
ments these studies may even result in the discovery of new phenomena which 
might give missing clues in which directions to search in the development of 
'new theories'. 

Because the content of the theory is essentially defined by Einstein's equa- 
tions it is not surprising that the main problem in the field is the analysis of 
these equations and their solutions. The following is meant to give an outline 
of the present situation and some insight into the questions to be dealt with. 
It is definitely not to be considered a survey of a field which has been growing 
since 1915. When I try to explain the situation I shall avoid technicalities to 
the extend to which this is possible without getting too vague in a field which 
abounds with technical questions by its very nature. References will be given 
mainly to illustrate a point or to direct the reader to more precise statements 
about concepts and results. The given references to specialized survey articles 
should allow the reader to get a more complete picture. 

2 Gravitational radiation, singularities, and black 
holes 

Einstein's analysis of gravitational radiation was based on the linearized field 
equations and the quadrupole formula. As indicated above, this gives essen- 
tially correct results in the 'weak' field of the binary pulsar and it is in fact still 
the basis for many calculations of gravitational wave emission in gravitational 
collapse scenarios. It can, however, hardly be expected to provide reliable an- 
swers in situations involving strong and highly dynamical fields, it is useless for 
calculating radiation generated by the coalescence of black holes. 

The reconsideration of the idea of gravitational radiation in the non-linear 
theory led in the 1960's to a concept of gravitational radiation, which does not 
require mathematical approximations but relies on the idealization of asymp- 
totic flatness in null directions f |24|). The latter is based on a picture of the 
overall behaviour of the gravitational field of an isolated self-gravitating system 
which assumes the field to become weaker and weaker in a characteristic way 
along null geodesies running out to null infinity. One might expect that 'close 
to null infinity linearized gravity takes over', but the situation turned out to 
be more subtle than that. Penrose has given a geometric and particular ele- 
gant characterization of asymptotic flatness in terms of the extendibility of the 
conformal structure through null infinity with a certain degree of smoothness 
f |l()5p . This suggestion associates with the far field of a selfgravitating isolated 
system a natural concept of radiation field but it poses at the same time difficult 
questions about the long time behaviour of gravitational fields and the nature 
of the field equations. It will be referred to in the following as Penrose proposal 
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(cf. injj, ISn] for further discussion and references). 

As it turned out, this was only a first reahzation of the importance of global 
considerations in the analysis of the non-linear field equations and the intro- 
duction of various general relativistic concepts. The dominating topics of the 
following years were concerned with physical situations such as gravitational 
collapse and cosmological singularities, were fields could be expected to show 
even unlimited growth over finite times. Some of the main results of that pe- 
riod of research were (i) the singularity theorems of Penrose |1(J6| and Hawking 
and Penrose , which showed that the occurrence of space-time singularities 
is a stable feature of gravitational fields, (ii) the analysis by Belinskii, Khalat- 
nikov, and Lifshitz j^, which led to the BKL conjecture about the general 
(oscillatory) behaviour of the field near cosmological singularities, and (iii) the 
development of the theory of black holes ([77], |H11)- 

It may be noted that 'singularity' is defined in the singularity theorems as 
the existence of a causal geodesic which is non-extendible and non-complete. It 
is left open, for which reason the curve should be non-extendible. In the case 
of specific explixit solutions more could be said concerning the nature of their 
singularities (cf. for the discussion of examples) but in general the methods 
given at the time were not sufficient to supply more information. 

Among the questions raised by gravitational collapse theory the most im- 
portant and still unsolved one is whether the evolution of gravitational fields 
admits a cosmic censorship principle ( |l(J7p which excludes in generic circum- 
stances the existence of naked singularities, i.e. singularities which could be seen 
by (possibly distant) observes. If such singularities existed stably under small 
perturbations this would reduce the predictive power of the theory. Giving a 
precise meaning to this principle is part of the problem. There exist different for- 
mulations in the literature such as the weak cosmic censorship principle, which 
asserts that singularities are hidden within black holes ((135|). and the strong 
cosmic censorship principle, which asserts that maximally extended space-times 
arising from generic non-singular initial data are globally hyperbolic (a notion 
explained below) ([102, [IMjl. 

There have been given fascinating accounts of these exciting developments 
by some of those involved (cf. [23, |31I, (SHli CSI]) and I shaU not try to 
repeat any of it here. Naturally, any research into a theory as complicated 
as general relativity starts from what could be called its 'fringes', defined by 
situations close to Newtonian ones, by static or stationary model situations, by 
configurations with other symmetry or other simplifying features, etc. While 
such studies, combined with perturbative calculations, have led to impressive 
results and far reaching extrapolations, the end of the story certainly still needs 
be told. Moreover, problems which may be considered by the early pioneers 
as having been settled a long time ago, may require new considerations in the 
darkness of new questions. 

The problem of cosmic censorship, questions about the strength of the sin- 
gularity, and the more specific but not unrelated questions raised by the BKL 
conjecture and the Penrose proposal represented the guiding problems for much 
of the subsequent research on the global evolution problem. 
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We may well only have scratched the surface of the domain of highly dynami- 
cal and strong field configurations and it is far from clear that we have exhausted 
the content of the theory. The general relativistic phenomena mentioned above 
had been predicted as theoretical consequences of the theory before they have 
been confirmed by carefully directed observations. Why should we have reached 
the end of it ? With our restricted theoretical and observational access to the 
domains alluded to above we may be missing out fundamental facts and unex- 
pected phenomena. The results by Choptuik ((SJ) on the phenomena occuring 
at the threshold of black hole formation may just give us a glimpse at things to 
come. 

The coalescence of black holes is expected to represent one of the strongest 
sources of gravitational radiation. Even in the clean cut pure vacuum situation 
still little can be said about this process under general assumptions, though 
reliable quantitative results are needed for analysing the data of the gravitational 
wave experiments. 

In recent years huge strides have been taken towards the goal of controlling 
the solutions to Einstein's field equations on large scale space-time domains. The 
results available now go, in any sense of the word, far beyond what was known 
around 1980. Nevertheless, getting qualitative and quantitative (resp. analytical, 
theoretical, and numerical) control on the long time evolution of gravitational 
fields under general assumptions is still the most important open problem in 
classical general relativity. 



3 The exploration of the solution manifold 

One of the main problems in controlling the behaviour of solutions to Eins tein' s 
field equations is posed by the Ricci-operator R^p [g] on the left hand side of Hl.l|) . 

In general coordinates (a;'^)A=o,...,3 the unknown metric g is given by a symmetric 
4x4 matrix (5^,y)^,i/=o,...,3 whose entries are functions of the coordinates. In 
terms of these unknowns the Ricci-operator reads 
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where we use the coefficients g^^ of the matrix {g^"^) inverse to (5^1^), which are 
of the form g^^ — {det{gf^^))~^ p^^' with polynomials p^^ of degree 3 in the g^^. 

In the case of a vacuum field with vanishing cosmological constant the matter 
fields and the energy-momentum tensor vanish and the various difficulties arising 
from matter models and matter equation arc not present. Equations Hl.l|l then 
simplifies to the vacuum field equations 



(3.2) 
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and the operator 1)3. l|l is all one needs to consider. But even somebody working 
in PDE theory for many years may just notice that H3.2|l is a equation of second 
order which is linear in the second derivatives of g^i, (i.e. quasi-linear), quadratic 
in the derivatives of and rational in the unknown g^^, but then may feel lost. 
Working out the content of the field equations without taking into account their 
geometric background is a hopeless task. Much of the richness of the system 
and important paths towards analysing the structure of its solutions will remain 
hidden. 

Usually one has to deal with in addition to equations Hl.lfl suitable matter 
models and the resultant matter equations. For physical reasons it is clearly 
most important to understand the behaviour of the resulting coupled systems. 
Matter equations import, however, their own specific difficulties (cf. IllOp . It 
is only because we do not want to be distracted by these individual properties 
that almost nothing will be said about results on solutions with matter fields. 
For several reasons the cosmological constant attracted increasing interest in 
recent years and there are available quite a few results on the corresponding 
solutions. Nevertheless, only the case A = will be considered in the following. 

Under general assumptions (no symmetries, algebraically non-restricted cur- 
vature tensors, no approximation requirements like low speeds etc.) qualitative 
results on the solutions to (|3.2II . in particular on large space-time domains, can 
only be obtained by applying geometric and abstract PDE methods to the analy- 
sis of the field equations, subject to various suitably chosen boundary conditions. 
Before formulating boundary value problems and proving anything, one needs 
ideas about fruitful problems and a hunch of what might be provable. Physical 
questions usually lead to deep mathematical problems whose analysis requires 
a considerable amount of ingenuity and the invention of new methods. But in 
the end they often allow for natural answers which illustrate the remarkable 
coherence of the theory. 

In the last 90 years there has been developed a host of methods to deduce 
information and heuristic results from the equations. The study of certain ex- 
plicit solutions exhibiting unexpected features, formal expansion type analyses 
of the solutions based on various representations of the field equations, studies of 
perturbations away from well- understood situations, and the use of topological, 
differential geometric, and PDE methods gave rise to far reaching generaliza- 
tions. 

In recent years numerical calculations have proven a powerful extension of 
the arsenal. A remarkable example is provided by the countable family of 
smooth, static, spherically symmetric solutions to the Einstein- Yang-Mills equa- 
tions discovered by numerical methods by Bartnik and McKinnon (jj^). These 
solutions certainly would not have been discovered by purely analytical methods 
up to this day (cf. (127^ for the complications of the analytic proof). 

A closer interaction between the analytical and numerical relativists appears 
to have a huge potential for further progress in the domain of relativity we are 
discussing here. In the following I shall try to point out some domains where 
such collaborations have been successful and also questions where a collabo- 
ration would be desirable and fruitful or even necessary. For a discussion of 
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questions more specific to numerical calculations I refer to the article and 
the literature given there. 

The abstract analysis of Einstein's field equations developed slowly. That 
these equations themselves satisfy the requirement of local causality was shown 
by Stellmacher (^128,) only in 1938. The first general local existence proof for 
solutions to Einstein's equations was given in 1952 by Choquet-Bruhat (ISni)- 
Until 1980 only existence locally in time was considered (cf. [31], jSHl for sur- 
veys). Choquet-Bruhat and Geroch obtained, however, an important uniqueness 
result, which also sheds light on the dangers of using the word 'global' ( 32 ). It 
states that an initial data set for Einstein's vacuum field equations determines 
a maximal, globally hyperbolic time evolution uniquely up to isometrics. 

Here a space-time (not necessarily solving any equations) (A/, g) is called 
globally hyperbolic if it contains a space- like hypersurface S such that each causal 
curve which cannot be extended as such meets S in precisely one point. In that 
case S is called a Cauchy hypersurface. Let (M, g) be a globally hyperbolic 
space-time and S a Cauchy hypersurface in M. Cauchy data for the linear wave 
equation Dg(j) = on then determine a solution which exists everywhere on 
M and which is unique. It will be shown below that the vacuum field equations 
can be understood as a system of (non-linear) wave equations. This makes the 
global uniqueness result plausible, but there are subtleties. 

There are known maximal globally hyperbolic solutions to Einstein's equa- 
tions which can locally be extended (as solutions) but the extension will not 
be unique. The hypersurface across which the extension takes place is called a 
Cauchy horizon. The presence of a Cauchy horizon signals a loss of uniquess, 
which is certainly something to worry about if the field equations are to predict 
the future. There arise delicate points here. So far we assumed everything to be 
smooth. Should we take the extension seriously if it can only be performed with 
low smoothness ? Which smoothness requirements will still make sense from the 
PDE point of view ? When do extensions still admit reasonable physical inter- 
pretations ? Where is the dividing line between 'extension of low smoothness' 
and 'weak singularity' ? 

One reason for the slow start of the general, abstract analysis of Einstein's 
equations is certainly the fact that already in the vacuum case the investigation 
of boundary problems for Einstein's field equations requires the study of four 
different differential systems: 

— the system of constraints, 

— the gauge system, 

— the main evolution system or reduced system, 

— the subsidiary system. 

The need for this is due to the diffeomorphism invariance of Einstein's field 
equations. In an initial value problem equation H3.2|) cannot have a unique 
solution in the sense of PDE theory. There exists a large set of diffeomorphisms 
which leave the metric invariant up to second order on the given initial time-slice. 
Therefore one has to take special measures to reduce a problem for Einstein's 
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equations to a standard PDE problem. Moreover, without further assumptions 
the operator (|3.1I) is not of the hyperboHc type which one might expect for an 
equation required to respect local causality. 

The systems above are not independent of each other. While some are fixed 
by the given boundary problem, there is a large freedom in the choices of others 
and changes of one of them ususally entails considerable changes in the others. 
A deeper understanding of their meaning, freedom, and interaction is required 
each time new problems are analysed. This will be illustrated below in the 
case of the subsidiary system. While it is of marginal interest in the analytical 
arguments, it seems to acquires an important role in the numerical construction 
of space-times. 

In the following the Cauchy problem will be considered most often. Of 
the many boundary problems studied for Einstein's equations, this is the most 
important but not the only relevant one. Under specific circumstances or as 
auxiliary problems other boundary problems may be equally important. To 
illustrate some of the questions which need to be considered in the context of 
existence theorems or in the numerical calculation of space-times, I begin with 
some remarks on the systems mentioned above and point out some old and new 
results. 

3.1 The constraint equations 

In the Cauchy problem for Einstein's field equations a prospective solution 
{M, g) of equation H3.2|l is characterized in terms of data prescribed on a 3- 
manifold S which is envisioned as being embedded as a space-like hypersurface 
into M. Since the initial data determine the solution near the initial hyper- 
surface uniquely, they contain the basic information on the solution. The con- 
struction, detailed understanding and interpretation of initial data is thus an 
important part of the Cauchy problem. Of the many results available now only 
a few can be discussed here and we refer to the recent survey |15| for more 
results and details of the methods only indicated in the following. 

It turns out that the data which need to be prescribed in a Cauchy problem 
for the vacuum field equations (|3.2|l are given by symmetric tensor fields hab 
and Xab which by the embedding of S into M will acquire the meaning of 
the Riemannian metric and the second fundamental form induced on S by the 
prospective solution g. As a consequence of the covariance of the field equations 
the data need to satisfy the vacuum constraint equations 

R[h]-XabX''' + {Xan'^0, D,Xa'-DaXc' = 0, (3.3) 

where D and R[h] denote the covariant derivative and the Ricci-scalar of the 
metric hab and the latter is used to move indices. These four quasi-linear equa- 
tions form an underdetermined elliptic system for the twelve components of hab 
and Xab- If one takes into account the freedom to perform transformations of 
the three coordinates and considers S as being determined in M essentially by 
the mean extrinsic curvature ip — ^ Xa"" , the rough function counting gives two 
degrees of freedom for the gravitational field. 
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The analysis of these equations depends on conditions which are not or only 
partially controlled by Einstein's equations such as the topology of S or the 
fall-off behaviour at infinities. Of particular importance are the cases where 
S is compact with arbitrary topology, where {S, hab) is asymptotically flat or 
asymptotically euclidean (corresponding to a hypersurface extending to space- 
like infinity where the metric hab approaches an euclidean metric), or where 
{S, hab) is hyperboloidal (corresponding to a hypersurface extending to null in- 
finity, like a space-like unit hyperbola on Minkowski space). In each of the 
non-compact cases one may consider fc > 1 asymptotic ends, so that for some 
compact subset if of S* the manifold S \ K has k components each of which is 
diffeomorphic to \ B where B is a. closed ball in R^. 

In these cases the construction of solutions to the constraints is well under- 
stood if the mean extrinsic curvature ip is assumed to be constant (-0 = in 
the asymptotically flat, tp = const. ^ in the hyperboloidal case). Following 
a suggestion by Lichnerowicz and using the behaviour of the equations under 
conformal rescalings of the metric, one finds that the metric can be chosen in 
the form hab = 'P'^ hab with some positive scalar function (j) which is to be deter- 
mined by solving some equation and a metric hab which is to be prescribed on S. 
If the symmetric trace-free tensor Xab then satisfies with respect to the metric 
hab the equation 5c Xa = the fields hab and Xab = 4>~'^Xab + "0 ^ab satisfy the 
vacuum constraints. This method, referred to as the conformal method, reduces 
the problem of solving the constraints to a linear elliptic system to obtain Xab 
and a decoupled semi-linear elliptic scalar equation for the conformal factor (j), 
called the Lichnerowicz equation. 

The solvability conditions for the latter in the compact and to some extent 
in the asymptotically flat case came along with the complete clarification of 
the Yamabe problem (a long standing mathematical problem, the final step of 
which was taken by Schoen who, remarkably, used ideas introduced by general 
relativity f |125| . cf. also ^95 )). The criterion is given in terms of the sign of 
the Yamabe number, an invariant of the conformal structure defined by {S, hab)- 
The general solvability condition in the asymptotically fiat case has been given 
only recently by Maxwell ([HHI)- It requires the (suitably generalized) Yamabe 
number to be positive. In the hyperboloidal case solvability conditions do not 
arise. 

The conformal method provides large classes of solutions to the constraint 
equations. In the case of non-constant mean extrinsic curvature "ip the conformal 
method does not lead, however, to a decoupling of the equations and the resultig 
simplifications. The existence of solutions to the constraints can still only been 
shown under severe conditions on ■0. The resulting restriction on the class of 
space-times which can be constructed from such data may be quite serious. 
Chrusciel, Isenberg, and Pollack have shown the existence of vacuum space- 
times which do not admit maximal (case -0 = 0) slices (03). Moreoever, time 
slices with ip ^ const, are encountered quite often in discussions of the evolution 
problem and the conformal method can not be applied to discuss the data 
induced on such slices. If such slicings are used in numerical studies it appears 
difficult to replace such data by improved data which are close to the given ones 
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and satisfy the constraints with higher aecuracy. 

We have restricted our discussion to the vacuum case for convenience only, 
there do exist methods to provide data for Einstein's equations coupled to var- 
ious matter fields. In those cases the data often have a direct physical interpre- 
tation. If the data comprise, for instance, a ball of perfect fluid with a vacuum 
exterior (cf. ||55| for a detailed discussion of this situation and its subtleties) 
we have a fairly clear idea about their meaning, though the large freedom to 
dispose of the exterior vacuum field still raises questions (cf. section 13. 5 .21) . 

In the pure vacuum case the physical meaning of the data is not so obvious 
and can in general hardly be assessed without analysing their evolution in time 
(cf. section [3.5. 2|l . Some interpretation if obtained for vacuum data with spe- 
cial properties. The singularity theorems and the cosmic censorship principle 
suggest that asymptotically flat data containing a trapped surface, an embed- 
ded surface S which is characterized by certain convergence properties of the 
out- and ingoing family of light rays orthogonal to E, develop into space-times 
containing event horizons and black holes. Evolving such data, which arise, for 
instance, if several asymptotic ends are present, is thus the usual method to 
model black holes and their coalescence. 

In numerical calculations, it may be advantageous to start from asymptot- 
ically flat initial data on a manifold S with an inner boundary S = dS which 
represents a trapped surface. Depending on the precise conditions to be achieved 
on S, the conformal method leads in this case to various overdetermined elliptic 
boundary value problems. Recently the nature of such problems has been anal- 
ysed and existence theorems have been proven by Dain (|51j). Dain, Jaramillo, 
Krishnan ([H]), and Maxwell ( 98 ). The initial data sets so obtained comprise 
exterior data which arise from non-trivial topologies as considered above as well 
as exterior data extending to data on S* = K'^ as described in section 13.5.21 

If asymptotically flat or hyperboloidal data are to be calculated numerically 
without imposing cut-offs at artificial finite boundaries, there arises the problem 
that the data can develop logarithmic sigularities at the asymptotic ends. In the 
asymptotically flat case conditions on the 'free data' under which logarithmic 
sigularities do or do not occur at space-like infinity have been discussed by Dain 
and Friedrich ([IS])- The analogous question for hyperboloidal data has been 
discussed by Andersson, Chrusciel, and Friedrich ([H]) and in great generahty 
by Andersson and Chrusciel (|S], [7]). 

Data with quite unexpected properties have been obtained recently by glu- 
ing techniques. A particularly remarkable idea has been introduced by Corvino 
in ,45,. The underlying method to exploit the underdeterminedness of the con- 
straint equations to obtain smooth, localized deformations of solutions to the 
constraints has been extended by Chrusciel and Delay f jifflj) and Corvino and 
Schoen ( I4(jj ). It provides a deeper understanding of the constraint map <&, 
which maps the fields hab, Xab onto the expressions on the left hand sides of 
equations H3.3|l and it allows one to construct solutions to the constraints which 
are not accessible by the conformal method. In particular, it enables one to 
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deform given asymptotically flat solutions to the vacuum constraints outside a 
given compact set to solutions which are exactly static or stationary in a neigh- 
bourhood of space-like infinity or to solutions which are asymptotically static or 
stationary at space-like infinity up to a given order or at all orders. 

As discussed below, the surprising freedom in modifying data in their asymp- 
totic domain sheds a new light on the Penrose proposal and it raises subtle con- 
ceptional questions about the calculation of wave forms characterizing isolated 
self-gravitating systems. 

Quite different aspects of initial data sets are addressed in the investigations 
of Penrose inequalities. With any asymptotically flat initial data set can be 
associated a certain invariant called the total mass or ADM-mass (^Sl)- It is 
obtained by performing a certain integral over a large sphere and taking its limit 
when this sphere is pushed to infinity so as to encompass the whole manifold S. 
In view of the weak conditions imposed on the initial data by the constraints 
it is quite remarkable that this mass could be shown to be non-negative and 
to be zero only for flat data. The Penrose inequalities may be considered as 
extensions of this positive mass theorem by Schoen and Yau and Witten 

the case of a space-like hypersurface S embedded in a space-time with 
event horizons so that the induced initial data set is asymptotically flat and S 
intersects the event horizon in a 2-surface E, these inequalities are expected to 
give a lower bound for the total mass of the initial data set in terms of the 
square root of the area of S. Remarkably, under the condition R[h] > on 
the Ricci scalar and certain assumptions which simplify the identification of the 
surface E, Penrose inequalities have been proven by Huisken and Ilmanen (|87|') 
and by Bray ( 26 ) (cf. |57] for a survey). 

A related class of problems is that of associating with an extended but finite 
space-time domain a notion of energy or energy-momentum. While some of the 
suggestions considered here played an important role in the above proofs of the 
Penrose inequalities, there does not exist a general agreement on the ultimative 
notions of quasi-local energy-momentum and other quasi-local quantities (cf. 
|129| for a detailed survey). 

It is interesting that Penrose arrives at the type of inequality named after him 
by invoking the 4-dimensional space-time picture and using a chain of arguments 
each of which raises questions itself f |l()8| 'l. He assumes in particular a version 
of weak cosmic censorship and makes use of the idea that after developing 
(something which is to become) an event horizon the space-time will settle down 
'in some appropriate but as yet ill-defined sense' to become a Kerr black hole. 
This idea is supported by the results on the black hole equilibrium problem (cf. 
j29| ) but a proof would require control on the long time evolution and estimates 
which describe in detail how the corresponding member of the Kerr family will 
be approached. 

Once the Penrose inequality can be derived by relying only on properties of 
initial data sets, one may ask whether the argument could be turned around 
and the inequalities or the techniques underlying their proofs could be used to 
obtain estimates to control the evolution of black holes. This will, however. 
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almost certainly require a proof of the appropriate Penrose inequalities under 
sufficiently general assumptions (including perhaps hyperboloidal data). As 
discussed above, general mean extrinsic curvatures create difficulties in solving 
the constraints. They also create difficulties in the present context. The first of 
equations H3.3|l shows that the condition R[h] > may be violated in the case 
of an unrestricted mean extrinsic curvature. New methods may be needed to 
obtain the inequality in such cases. 

3.2 The gauge system 

To apply PDE techniques to the local evolution problem, one has to impose 
gauge conditions, restrictions on the freedom to perform diffeomorphisms or co- 
ordinate transformations. In the numerical calculation of space-times a number 
of unsolved questions are related to the gauge problem. To illustrate the general 
argument and related problems without being too vague, 1 indicate one specific 
reduction procedure by which the initial value problem for Einstein's field equa- 
tions is cast into a Cauchy problem for a hyperbolic system. The one chosen 
here yields the most concise expressions. 

In the given local coordinates x'^ on M the expression ^'^.1^ can be rewritten 
in the form 

R^.. = 5^" g>..,xp + V(^r,) + Ta V g^s g^" T/, + 2Ts\g'P g^^ T,) " p. 

(3.4) 

Here the comma indicates partial derivatives, the V,j^n = ^/"^ g^^ {gx-n.v + 
gvX.ri — gvq.x) are the Christoffel symbols, V is the covariant derivative operator 
of (7^^, and the summation rule applies. The F'' denote the contracted Christoffel 
symbols = g""* F^ ^ which (together with the functions F^, = g,^^ consid- 
ered in the following) are being formally treated as if they defined a vector field 
(which, of course, they do not). Thus F^ = g^f^ F^ and V^F^ — df_,T^—Tf^ ^ ^ V\. 

The form %\A\ emphasizes the first term on the right hand side of (|3.1|l 
which is obtained by applying to the unknown g^i, a wave operator, a type of 
differential operator for which a good theory is available. The following three 
terms of second order in (|3.1|l prevent the direct application of PDE results. 
They are hidden in the representation (|3.4() in the second term on the right 
hand side. It turns out that the apparent difficulties dissolve once the role of 
the contracted Christofffel symbols F"^ is recognized ( 63 ). 

Let {M,g) denote some Lorentz manifold and let S = {t = 0}, with some 
coordinate function t, be some space-like hypersurface of it. Consider a map 

9 x'^ ^ F'-'-{x^) E Kf^. Ignoring subtleties arising from differentiability 
questions, everything is assumed to be smooth. Local Cauchy data x^, dtx^ on 
S determine a unique local solution to the Cauchy problem for the semi-linear 
system of wave equations 

Ugx'' = -F^'{x^), (3.5) 

where Dg — V denotes the scalar wave operator defined by g. If the dx^ 
are chosen initially to be pointwise linearly independent, the solution provides 
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a local coordinate system x'^. In terms of these coordinates the system above 
simply takes the form 

r^(a;^) = F^(a;^). (3.6) 

As a consequence 

— by a suitable choice of coordinates the contracted Christoffel symbols 
can locally be made to agree with any prescribed set of functions F^^ , 

— in turn, these gauge source functions and the initial data determine 
the coordinates uniquely, 

— for a given metric g any coordinate system is characterized by suitable 
gauge source functions and initial data. 

This suggests replacing the functions in H3.4|) by freely chosen gauge 
source functions F'^. The vacuum field equations then take the form 

= Rf,, = -^9^' .9m-.Ap+V(^ F,) " ^ g,^s g^' T/ , + 2Ts\ ^^(^ T,) " p, 

(3.7) 

of a system of quasi-linear wave equations for the g^,j, which represents the 
main evolution or reduced system of our procedure. For this system the local 
Cauchy problem with appropriate data on a space-like hypersurface S is well 
posed, which means that there can be shown the existence and the uniqueness of 
solutions and their stable dependence on the initial data (cf. [75] for a detailed 
discussion). A few interesting observations to be made about this procedure. 

The metric coefficients obtained as solution to (|3.7II are given in terms of 
coordinates x^^ which are determined implicitly by the gauge source functions, 
the initial data, and equation H3.7|) . Do these coordinates really satisfy the gauge 
system ()3.5|l resp. its implicit coordinate expression (|3.6|l ? For the moment we 
assume this to be the case and consider this question again in section (|3.4f) . 

The domain on which the x'^ form a good coordinate system depends on 
the initial data, the gauge source function, and on the solution g itself. Since 
information on g is only acquired by solving H3.7|l . little can be said a priori on 
the domain of existence of the coordinates. In this respect there is in general no 
difference between harmonic coordinates (more appropriately called now wave 
coordinates), characterized by — 0, and coordinates defined by other gauge 
source functions. Without special precautions nothing prevents the dx^ to be- 
come linearly dependent, the slices x^ = t = const, with {t = 0} = S" to turn 
time-like, the coordinates to develop an undesirable asymptotic behaviour, etc. 
If a coordinate system turns bad, one will have to construct further coordinates 
and may end up with a collection of overlapping coordinate patches which serve 
to define the manifold M underlying the solution space-time. 

In practice, in particular in numerical calculations of space-times, one would 
like to avoid such situations and try to find coordinates which cover the entire 
solution. An interesting way of using gauge source functions to control the 
evolution of the slicing in a numerical code has recently been put forward by 
Pretorius in his work^ following up |112| . It remains to be seen whether the 

^In preparation 
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method will allow him to steer the slicing unscathed through all the dangers of 
a long time evolution. 

Lindblad and Rodnianski discuss the small data, global existence for the 
Einstein vacuum equations in wave coordinates and it turns out that the coor- 
dinates can be arranged to cover the complete solution space-time ((ST])- These 
solutions have an asymptotic structure qualitatively similar to that of Minkowski 
space. If the data were slightly changed one would expect the coordinates to 
be still well behaved by stability considerations. Such considerations become 
quite delicate in global problems, however, and unexpected things may happen 
when the data are increased to admit the development of singularities and black 
holes. In any case the behaviour of the coordinates has to be controlled in the 
context of the evolution problem, in which little is known about the metric a 
priori. 

In [2] have been discussed coordinates which cover the entire Schwarschild- 
Kruskal space-time up to the singularity and even extend smoothly through 
the null infinities. Again, it is unclear how they will behave if the underlying 
solution space-time is perturbed. This general class of coordinates is based 
on certain geometrically distinguished curves called conformal geodesies and 
the defining equations of these curves contain elements which might allow one 
to control the behaviour of the coordinates over long time intervals. While 
these coordinates can be characterized in principle in terms of the gauge source 
function considered above, there is no way of identifying these functions without 
knowing the solution g. Therefore one needs a different reduction procedure to 
incorporate these coordinates into hyperbohc evolution equations In 
this reduction the coordinates are characterized by explicit conditions on the 
unknowns and one might hope that the question which we left open above does 
not arise. It comes back as the question of certain constraints being satisfied 
during the evolution. 

The last example shows that it is not easy to identify good gauge source func- 
tions. It shows also that there may be good reasons for analysing reductions 
different from the one indicated above. The role of the gauge source functions 
may be assumed then by quite different quantities. Motivated by problems in 
numerical relativity there has been considered a large variety of different reduc- 
tions in recent years, based on different representations of the field equations, 
different unknowns, and different types of gauge conditions. Many of these 
reductions lead to hyperbolic main evolution systems. 

3.3 The main evolution system 

The main evolution system is clearly most important for working out local 
existence, uniqueness, smoothness, and more specific properties of solutions and 
for calculating solutions numerically. To avoid entering technicalities I shall only 
make a few general remarks about it. 

In analytical work the main goal in choosing this system is to be able to 
exploit the intrinsic hyperbolicity of the equations. This does not mean that 
the system needs to be hyperbolic. Some useful gauge conditions are elliptic 
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in nature and there have been studied for instance hyperbolic-elUptic reduced 
systems by Andersson and Moncrief CfTT]'). 

In recent years there has been a tendency in the general investigation of 
non-hnear evolution equations to study solutions of low smoothness ( |13Up . 
Klainerman and Rodnianski (|93j) study Einstein evolution under smoothness 
assumptions which are weaker than those considered up to a few years ago (cf. 
|SBj) and Maxwell (HU) discusses the existence of rough solutions to the con- 
straints equations. If one is mainly interested in physical phenomena, these 
activities may appear quite esoteric. As indicated above, questions about pre- 
cise and low smoothness requirements can become, however, inavoidable in the 
discussion of singularities and Cauchy horizons. Moreover, trying to push to 
their lower limits the smoothness requirements under which the Einstein equa- 
tions still make sense forces one to explore the specific structure of the equations 
much more carefully and it is bound to lead to more precise information on the 
evolution. 

Once local existence has been treated, there may be used other methods to 
get control on the long time evolution. While some hyperbolic main evolution 
systems imply energy estimates involving the Bel-Robinson tensor, this tensor 
may be used independent of any reduction to derive estimates on the solution. 

In any case there needs to be erected some kind of rigid space-time structure, 
a foliation by space-like or null slices or a fixed coordinate system, relative to 
which estimates of the metric field are expressed and the evolution of this struc- 
ture itself can be controlled. In the case of the global or semi-global non-linear 
stability results mentioned below one can lean back on information supplied 
by an explicit reference solution, choose on it a foliation determined by some 
suitable (evolution) law, and construct a foliation on the perturbed solution gov- 
erned by the appropriately perturbed law. If the coalescence of two black holes 
is to be modelled under general assumptions, however, there are no reference 
solutions available and one has to develop an intuition for foliations with long 
life times and good evolution properties in the context of an existence proof or 
in the course of an numerical calculation. 

These remarks also show that the needs of analytical and numerical work are 
different. In the latter one has to rely on an explicit main system for all times of 
the evolution (or at least change the system only a finite number of times). For 
the choice of this system the main requirement is the stability of the numerical 
evolution. So far this has essentially been a matter of trial and error and there 
appears to be no way to translate this requirement into a precise criterion in 
terms of algebraic or other properties of the system. Manifest hyperbolicity is 
reasonable but apparently not sufficient for that purpose. There is in use a class 
of systems, the so-called BSSN systems (QHIj cf. also P, [7S], fTS] for hyperbohc 
versions) some of which are not manifestly hyperbolic but nevertheless seem to 
lead to stable numerical evolutions. Understanding why this should be so or for 
which class of problems stability fails for these systems is a theoretical challenge. 
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3.4 The subsidiary system 

In our discussion of the gauge system we left open the question whether the 
imphcit gauge condition is preserved under the evolution defined by H3.7|) . 
Here is the analytical argument. Equation (|3.7|l is of the form 

-R^ii/ = Qi,), (3.8) 

where = T — with the T^^ calculated from the solution g^i,. The require- 
ment that — may be considerd as a kind of constraint for H3.8|l (it is in 
fact related to the constraints considered in section (|3.1|l '). The twice contracted 
Bianchi identity, which holds for any metric, reads V^(i?p,y — ^ i?.g^,y) = 0. Ap- 
plying this to the equation above gives 

+ i?^ . g,, = 0, (3.9) 

and thus a subsidiary system of wave equations for the quantities Q^. 

It turns out that for data satisfying the constraints and the gauge condition 
= on the solution to H3.7|l satisfies 9tQ^ = on S*. Since we have 
thus vanishing Cauchy data for the hyperbolic system H3.9|l . it follows that 
vanishes and g^j/ does solve (|3.2|l . This closes the argument in the continuum 
model. Note that only the homogeneity and the uniqueness property of (|3.9|l 
are being used here. There is no further role for (|3.9(l in analytical studies. 

The situation is quite different in the discrete model. Most numerical calcu- 
lations of solutions to Einstein's equations are being plagued by an undesirably 
fast growth of constraint violations. In fact, many workers in the field report 
on seemingly unmotivated catastrophic blow-ups of numerical calculations at a 
stage of the numerical time evolution where coordinate or curvature singulari- 
ties were not to be expected. Understanding this situation requires and deserves 
a major effort. In the following I shall not discuss any of the remedies which 
have been suggested (cf. [SH]? IHODj or the stabihty analyses of the subsidiary 
system, considered as a linear system on a given background (cf. jHll, 1138) 1. I 
would like to identify instead possible sources of the problem in the analytical 
structure of the equations. 

Analytically, 1)3.9)1 is simply a differential identity implied by ()3.8)l . It is hard 
to see, however, how to devise a numerical scheme for the second order wave 
equations 1)3.8)) for which the subsidiary system 1)3. 9)1 . which is of third order in 
the metric, could be identified as an identity. The relations between the two 
systems will therefore become obscured, and the development of the constraint 
violation is not easy to analyse. 

If we observe ()3.8ll in ()3.9)l , the latter takes the form 

V^V^Q, + g^V(^Q,) =0, (3.10) 

of a manifestly non-linear wave equation. In the continuum model this equation 
would still imply = 0. In the discrete model the quantity comes, how- 
ever, with an error initially and develops further errors during the evolution. 
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The detailed propagation properties of (|3.1U|I therefore become iraportant. If 
H3.9|l were a Unear system on a given background, standard energy estimates 
would admit and cannot exclude an exponential growth of the unknown Q^i but 
would admit nothing faster than that. The non-linearity of (|3.1U|) might induce, 
however, a much faster growth of the constraint violation. 

Analysing this situation is not easy. While H3.8|l can be studied indepen- 
dently, equation (|3.10() does not decouple from H3.8|l . The latter supplies the 
metric defining the background for H3.10|l and the growth of has an effect 
on the evolution of the metric by (|3.8|) . The coupling of (|3.1U|) to (|3.8|l defines 
a second, though less direct, non-linearity in the evolution of Q^. Neverthe- 
less, considering the background metric in (|3.1U|) as given, should provide some 
understanding of the propagation properties of that equation. 

If one considers (|3.10l) as an equation on Minkowski space, it turns out that 
it is easy to find Cauchy data and dtQ^ on {t = 0} for which the solutions 
develop poles after finite coordinate times > (cf. [JH] for details). These 
data can be chosen pointwise as small as one likes, though ^ oo if the data 
approach zero. This suggests that one may have a relatively stable numerical 
evolution for a while but at a certain stage effects due to the non-linearity in 
(|3.10|l take over and induce a catastrophic collapse of the calculation after some 
finite time. 

There are quite a few interesting questions to be answered and the situation 
certainly deserves closer analytical and numerical study. We have discussed 
here the non-linearity of the subsidiary system in a specific example. Most 
likely, non-linearities are found in all subsidiary systems. To what extent their 
effects are different and whether there exist preferred cases remains to be seen. 

3.5 Global studies 

In the following we want to present some global or semi-global results on the 
existence of solutions to Einstein's field equations. The field has seen a rapid 
development in recent years. It will not be possible here to give due reference to 
all the important contributions and we will have to make a somewhat arbitrary 
choice. In the following it will be convenient to make a distinction between 
vacuum solutions which are asymptotically flat and solutions with compact time 
sections. In spite of the fact that there are important cosmological models 
with non-compact time slices, by a cosmological space-time will be meant in the 
following always a vacuum solution with compact time slices. 

3.5.1 Cosmological space-times 

In the context of cosmological solutions there exists a large variety of possible 
assumptions on the (local) symmetries and the topologies of space-sections. The 
tendency of much of the recent work has been to analyse the global behaviour 
of the solutions under strongly simplifying assumptions and then, relaxing them 
step by step, to work ones way up to develop in the process the insight and the 
technical means to analyse in the end also quite general classes of solutions. Of 
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the many results obtained in this program only a few can be considered in the 
following. 

What are the questions to be asked ? If one tries to work out a 'global' 
existence result by analysing a Cauchy problem for a certain class of solutions 
and it can be shown that the maximal globally hyperbolic solutions determined 
by the data are all geodesically complete in the past and in the future one 
may still be interested in the precise asymptotic behaviour near past and future 
time-like infinity but the essential goal has been reached. 

In general the situation will be more complicated. It may happen instead 
that the solution will be geodesically complete in one time direction (or not at 
all) but in the other direction there are obstruction to completeness. Extreme 
(and, in a sense, the clearest) possibilities which can occur are that the solu- 
tion approaches a smooth Cauchy horizon or that it develops of a curvature 
singularity at which an invariant built from the curvature tensor such as the 
Kretschmann scalar R^t^Xp R'^"^'' becomes unbounded. This may happen glob- 
ally when a family of space-like hypersurfaces approaches the corresponding end 
of the evolution or it may happen only locally. There are possible all kinds of sit- 
uations in between and combinations thereof. One would like to have a detailed 
description of the respective behaviour and relate specific types of behaviour to 
the 'size' of the corresponding subset of the space of initial data. 

Since any statement about cosmic censorship must refer to generic classes of 
initial data, the analysis of solutions with symmetries cannot give a final answer 
to the question of cosmic censorship. It can, however, provide insight into the 
nature of the problem to analyse the question in the class of space-times with a 
given symmetry. Moreover, if the subset of solutions which do not admit Cauchy 
horizons contains a dense open subset of the data set one would be prepared to 
consider this as evidence for cosmic censorship (restricted cosmic censorship). 

The 'simplest' and still quite interesting case to consider is given by space- 
times which are spatially homogeneous in the sense that they admit a 3-dimensional 
space of Killing fields which at each point generate the tangent space of a fo- 
liation by (compact) space-like hypersurfaces. In the universal covering space 
the (not necessarily compact) leaves of the foliation can then be considered as 
orbits of a 3-dimensional Lie group G of isometries. These Bianchi space-times 
have been classified in terms of the Lie algebras of their isometry groups. 

For spatially homogeneous space-times the Einstein vacuum equations es- 
sentially reduce to systems of ODE's. Nevertheless, their analysis turns out 
quite difficult, not because of the various cases in the Bianchi classification but 
because of the various phenomena which can occur. The Bianchi classifica- 
tion consists of two classes. We shall consider only one of them, in which the 
structure constants of the Lie algebras can be specified in terms of a symmetric 
matrix which can be assumed to be diagonalized (case A) . In that case the Ein- 
stein vacuum equations have been written by Hsu and Wainwright f |134| ') as an 
ODE and a constraint for an unknown u = (Ni, N2, N^, E+, E_) of the form 

= f{u), q{u) = 1, 
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with q a quadratic polynomial in the components of u. The first three compo- 
nents of u contain essentially the information on the Lie algebra and the other 
two information on the second fundamental form on the leaves of the foliation. 

When the Lie algebra is cummutative (Bianchi I) only E_ do not vanish 
and the constraint reduces to the equations + = 1, which defines the 
Kasner circle parametrizing these solutions. The fields E+, S_ are in fact 
constant, i.e. the solutions are fix points of the dynamical system above. The 
solutions are given explicitly by the Kasner metrics 

a=l 

on ]0, oo[xM^ or, after a periodic identification, on ]0, oo[xT^. The real numbers 
Pa satisfy the Kasner relations X]i=iPa ~ ^ ^^^'-^ "12^=1 Pa — 1 ^^'^ related to 
the two unknowns above by E_|_ = |(p2 + P3) — 1 and E_ — ^{p2 — Pa)- 

The three points (—1,0), (^,±-^) on the Kasner circle correspond to flat 
solutions, which in the case without identifications are isometric to open subsets 
of Minkowski space given by the future of the intersection of two null hyper- 
planes. In the non-flat case the causal geodesies are incomplete in the past 
(where t I 0) and approach a curvature singularity there. Note that the volume 
of the time slices grows monotonically with t but the growth happens locally in 
an anisotropic way because the Pa can not all have the same sign. 

Of particular interest is the case where the group G is SU{2) (Bianchi IX). 
It contains a subclass, the Taub type IX solutions, characterized by linear con- 
ditions on the components of u. These solutions can also be given explicitly 
and it turns out that they admit in the past (compact) Cauchy horizons and 
smooth extensions (there do exist in fact non- isometric extensions ( .40 )) across 
it, resulting in the Taub-NUT space-times. 

The remaining 'generic Bianchi IX solutions' or proper Mixmaster solutions 
( \1{)V\) show an interesting behaviour. They have a curvature singularity which 
is approached by the solutions with an oscillatory behaviour. The projections 
of their trajectories into the (S+, I]_) - plane moves into the Kasner circle and 
stays there, approaching subsequently different points on the circle (cf. 22 and 
|115j . |122| for numerical and analytical studies). The detailed studies of the 
Bianchi models show that curvature blow-up and the non-existence of Cauchy 
horizons in the contracting direction is a feature of generic Bianchi space-times, 
which shows that retricted strong cosmic censorship holds in this class of models 
([ng, cf. also 03). 

While being interesting in themself, the results on the Mixmaster solutions 
are expected to be of a much wider significance. The BKL conjecture mentioned 
in section El suggests that the behaviour of these solutions near the singularity 
provides a model for the local behaviour of general cosmological solutions in the 
neighbourhood of singularities. 

Moncrief initiated with the article |101) the global study of a class of vac- 
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uum solutions with two commuting space-like Killing fields and compact space 
sections called Gowdy space-times. With the spatial topology being that of the 
3-torus the metric can be written in the form 

g ^ Vre^ {-dt^ + dx^) + t {e^{dy + Q dzf + e^^ dz^), 

with X, y, z each being a coordinate on the circle 5^, the coordinate t taking 
values in M-|_, and the coefficients A, P, Q depending only on t and x. 
For these metrics the Einstein vacuum equations read 

Pm - Pr. - -\ P.t + e^iQ^t - Q%), 

Q ,tt ~ Q,xx = ~— Q,t — 2 {P^t Q,t — P,x Q,x): 

and 

X,x = -2t{P^tP,x + e^''Q,tQ,x), 

If Cauchy data for P and Q are given for which the integral over the circle 
of the expression on the right hand side of the third equation vanishes, the 
discussion reduces to the analysis of the semi-linear system of wave equations 
for the functions P and Q, which essentially represent the two polarizations 
states of the gravitational field. 

After Moncrief had given the first global existence proof for the solutions, 
these solutions and the nature of their singularities were studied by many au- 
thors. The analysis was largely assisted by the numerical work initiated by 
Berger and Moncrief ([ini). It showed that the solutions tend to develop strong 
gradients (spikes) in the approach towards the singularity, which brought an 
important aspect into the analytical discussion. Kichenassamy and Rendall 
constructed families of solutions with singularities by Fuchsian methods (cf. 
the survey |119| and also the work by Chae and Chrusciel (Eni))- Rendall 
and Weaver constructed families of Gowdy solutions with spikes from solutions 
without spikes, and distinguished 'true' spikes, which have a geometric meaning, 
from 'false' spikes f |12Up . Building on this and earlier work, Ringstrom recently 
showed that for a 'generic' set of initial data the corresponding solutions exhibit 
a curvature blow up on dense open subsets of the singularity f |128j V Combin- 
ing this result with the work by Chrusciel and Lake on the occurrence Cauchy 
horizons in Gowdy space-times he concludes that this class of space-times 

satisfies a restricted strong cosmic censorship principle. 

This short discussion hardly gives credit to the many important contributions 
which led to the 'final' answers. A more complete picture of these interesting 
developments can be obtained from the articles 0, iS^, and ,21^. The latter is 
also particularly interesting because it highlights the remarkably successful and 
still ongoing interplay between numerical and analytical studies in this field. 

The results mentioned above do not finish the analysis of solutions with 
two Killing fields. In fact, the Gowdy metrics considered here only define a 
'negligible' subset of the set of all solutions with two Killing fields. 
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As a further step in the program solutions with only one Killing fields should 
be analysed. A semi-global, non-linear stability result in this directions is ob- 
tained by Choquet-Bruhat, who generalizes previous work with Moncrief to 
show future completeness for a class of i7(l)-symmetric vacuum solutions on 
manifolds of the form M = R x E x 5^ with Cauchy hypersurface diffeomorphic 
to E X S*^, where E is an orientable, compact surface of genus greater than 1 
and the space-like Killing fields are assumed to be tangent to the fibres of the 
fibration defined by the projection M — > R x E (1311). 

Generalizing even further, Andersson and Moncrief obtain a semi-global, 
non-linear stability result for vacuum solutions without imposing any symmetry 
conditions (^2) (the solutions of |3T] are not included). Denote by V the 
interior of the future light cone at the origin in Minkowski space and by r 
the Minkowskian distance from the origin. Identification of points of V by the 
action of a suitable discrete subgroup of the Lorentz group yields a reference 
space-time of the form 

M R X 5, g = -dr^ + h, 

where (S*, h) denotes a 3-dimensional, compact hyperbolic space of sectional 
curvature —1. We assume dr to be future directed. The authors consider cases 
where {S, h) satisfies a certain rigidity condition. Identifying 5* with the set 
{r = 1}, the metric above induces on S the reference data {hab,Xab = hah)- 
It is shown in |12j that vacuum data on S sufficiently close to rigid reference 
data develop into solutions of the vacuum field equations for which the causal 
geodesies are future complete. The authors also desribe the asymptotic decay 
of their solutions towards the reference solution. 

For recent attempts to control the behaviour near the singularity under 
'general' assumptions, to give precise meaning to the BKL conjecture, and to 
develop tools which would allow one to decide on its valitity, we refer to the work 
by Andersson et al. ( 10 ) and Garfinkle ([ZEl)- There appears to be a general 
expectation that the BKL conjecture will turn out to be basically correct. 

3.5.2 Asymptotically flat space-times 

Asymptotically flat space-times provide the basic model of isolated self-gravita- 
ting systems such as stars, star systems, black holes, etc. and as such they 
are important for discussing many observable general relativistic phenomena, 
the analysis of radiative phenomena being at present the most important and 
urgent one. The modelling of stars is clearly an important task. Nevertheless we 
shall concentrate again on the vacuum case or situations with field theoretical 
matter models. On the one hand, progress in these cases will provide important 
insights into the behaviour of dynamical black holes, the coalescence of black 
holes, the resulting radiation fields etc., and possibly quantitative results about 
the latter, while on the other hand even these configurations still present major 
challenges. 

Asymptotically flat space-times pose problems which do not occur in the 
cosmological context. The space-like slices are of infinite extent and on an 
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asymptotically flat slice the fall-off behaviour near space-like inflnity can, in 
terms of coordinates which realize the asymptotic flatness conditions, not be 
faster then 

772 1 

5^1^ - '^P'' = ^-|^ (^M"^ + as ^ oo, 

unless the total mass 777 vanishes, in which case the space-time is flat by the 
positive mass theorem. This poses major problems for global existence proofs 
and, in particular, for the verification of the Penrose proposal under general 
assumptions. 

The infinite extent of space-like slices also requires additional considerations 
in the numerical study of space-times, because numerical calculations need to be 
performed on finite computational grids. Quite a few technical considerations 
and also conceptional questions, some of which will be indicated below, are 
related to that fact. 

Since gravitational radiation can either escape to null infinity or fall into 
a black hole, which is possibly (cf. the discussion below) generated by the 
radiation itself, the completeness properties of the space-times in the future 
and the structure of their time-like infinities can vary considerably. Already in 
stationary examples such as the Schwarzschild-Kruskal, the Reissner-Norstrom, 
and the Kerr solution there arise extreme and quite distinct situations near 
time-like infinity, where singularities, event horizons, null infinities, and Cauchy 
horizons seem to meet in the standard causal pictures (|H21)- Hardly anything 
is known about the possibilities under general assumptions. 

The central question in the field is again whether cosmic censorship is a 
valid principle. Since any precise statement about that principle will use the 
word generic, the general answer can only be obtained by analytical methods. 
Numerical methods are likely to play, however, a major role in the investigation 
of the possibilities. 

At present the main question related to physical observations is concerned 
with concepts of radiation, the precise asymptotic behaviour at null infinity, 
and methods to derive quantitative results about radiative properties. Here 
numerical methods are bound to play a dominant role in the end, but before 
that quite a few analytical questions will need to be answered. 

The two questions above are not independent of each other. Even if one 
only seeks to calculate the radiation field, if one were not interested in the inte- 
rior of black holes and the structure of singularities, and if one knew somehow 
that singularities were always hidden behind event horizons, one could hardly 
ignore in the analysis of the long time evolution of gravitational fields the ten- 
dency of solutions to develop singularities. Furthermore, since the location of 
an event horizon is not known at a finite stage of the evolution, the domain of 
outer communication, comprising the far fields, can in the usual approach to 
initial value problems not be cleanly separated from the interior of a black hole. 
Present attempts in numerical calculations to cut out the singularity from the 
computational domain may become a delicate matter in long time calculations. 
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In the case of asymptotically flat solutions the number of interesting simpli- 
fying assumptions is much smaller than in the cosmological context. If station- 
ary solutions are excluded and one insists on complete and regular far fields, 
only spherical or axial symmetry remain as admissible symmetries. Additional 
fields need to be considered to analyse any dynamical behaviour in a spherically 
symmetric setting, because spherically symmetric vacuum solutions are, by the 
Birkhoff theorem, locally isometric to patches of the Schwarzschild-Kruskal so- 
lution. Of the large body of work dealing with spherically symmetric situations 
the following three contributions are particularly important. 

In a remarkable series of articles Christodoulou analyses the formation of 
black holes and singularities for the spherically symmetric Einstein-scalar field 
equations. He gives conditions on the data for the avoidance and for the devel- 
opment of singularities respectively. He further gives conditions under which the 
singularity will be hidden by an event horizon but he also finds solutions whose 
singularities can be seen by distant observers, thus showing the occcurence of 
naked singularities. Finally, he shows that the existence of naked singularities is 
in fact an unstable property for the spherically symmetric Einstein-scalar field 
equations, which supports the cosmic censorship hypothesis (cf. [35] and the 
references given there). 

In a seminal article (1211) and subsequent work Choptuik studies numerically 
one-parameter families of spherically symmetric Einstein-scalar fields for which 
the solutions disperse for values of the parameter below but form black holes for 
values of the parameter beyond a certain threshold value. This allows him to 
discover phenomena which would have been difficult to find by purely analytic 
methods. In particular, he finds self-similar critical solutions, the vanishing of 
the black hole mass as the parameter approaches its critical value, and a certain 
scaling of the black hole masses with universal critical exponents. This work 
initiated quite a number of further investigations involving various different 
matter models and different types of non- linear equations (cf. [7^1 for a survey). 

Dafermos and Rodnianski study the spherically symmetric Einstein-Maxwell 
scalar field equations, assuming that a regular event horizon has formed ('|49|'). 
They give rigorous proof to Price's result f !113j) that perturbations of gravi- 
tational fields show in terms of a suitably chosen advanced time coordinate a 
polynomial decay on the event horizon near time-like infinity. Moreover, they 
confirm results by Israel and Poisson concerning the occurrence of weakly 

singular Cauchy horizons (cf. also |l()3j for the general background and 
PT] for numerical studies of the interior of black holes). The significance of 
this result concerning the question of strong cosmic censorship remains an open 
question as long as only spherically symmetric situations are considered (fortu- 
nately some kind of censorship is enforced already by the presence of the event 
horizon). 

Compared with the richness of these detailed results and observed phenom- 
ena, the study of the large scale structure of asymptotically flat vacuum solutions 
without symmetries is still in its infancy. Compared, however, with the technical 
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difficulties to be overcome, quite some progress has been achieved in the last 
twenty years. All the global or semi-global results concerning asymptotically 
flat solutions without symmetries which have been obtained so far are stability 
results which show the existence of geodesically complete or future complete 
solutions for small data, i.e. for data which are in suitable norms close to the 
data of a well underderstood reference solution, provided by Minkowski space 
or parts of it. Analogous results for comparison solutions with black holes are 
not available yet. 

An early semi-global existence result was obtained in the article where 
it was shown that smooth hyperboloidal data sufficiently close to Minkowskian 
hyperboloidal data develop into solutions of the vacuum field equations which 
have a smooth future complete structure at null infinity for which time-like 
infinity is represented in suitable conformal extensions by a regular point. The 
analysis uses the conformal behaviour of the Einstein equations in an explicit 
way and it is carried out in terms of a conformally rescaled metric with respect 
to which null infinity and time-like infinity are at a finite location. The solutions 
show in particular the peeling behaviour, so that along outgoing null geodesies 
with affine parameter r the components of the conformal Weyl tensor satisfy in 
a suitable orthonormal frame a fall-off behaviour near null infinity which can be 
expressed in terms of certain entire powers of r^^. 

The particular initial value problem was considered to avoid in a first step 
the complications at space-like infinity. While the setting was thus intended as 
a preparation for a global study, it acquired in the meantime some interest for 
the numerical calculation of radiation fields based on the underlying conformal 
field equations (ED, lEHI, IHHI)- 

Soon afterwards attempts were made to use this result to construct complete 
solutions. The idea was to avoid the complications at space-like infinity by 
constructing Cauchy data for which the evolution near space-like infinity could 
be controlled explicitly and the existence of smooth hyperboloidal slices close to 
Minkowskian ones could be shown. Cutler and Wald managed to construct non- 
trivial data for the Einstein-Maxwell equations which are spherically symmetric 
outside a compact set (02]). Thus they were able to establish the existence of 
geodesically complete solutions to the Einstein-Maxwell equations with smooth 
and complete structures and non-trivial radiation fields at future and past null 
infinity (cf. More recently, using the result of Corvino mentioned above, 

Chrusciel and Delay showed the existence of a large class of vacuum solutions 
with these properties ([3H|)- To some extent this justifies the Penrose proposal 
but the fact that these solution are exactly Schwarzschild near space- like infinity 
clearly leaves space for generalizations. 

The first global existence result for a general class of asymptotically fiat 
data was obtained by Christodoulou and Klainerman (IMI)- They exploit the 
conformal behaviour of the solutions only indirectly. With a considerable effort 
they skillfully manage to control the behaviour of the fields near space-like in- 
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finity. They show that vacuum data satisfying certain regularity conditions near 
space-hke infinity (which include the vanishing of the linear momentum) and 
suitable smallness conditions on the initial slice develop into a unique, globally 
hyperbolic solution of Einstein's vacuum field equations which is geodesically 
complete and asymptotically flat in the sense that the Riemann curvature ap- 
proaches zero along any geodesic if the afline parameter tends to infinity. The 
global structure of their solutions is qualitatively similar to that of the space- 
times considered above, with no further conditions on the data, however, the 
peeling behaviour in its usual form is not satisfied and the smoothness of the 
conformal structure at null infinity is thus weaker than that required by the 
Penrose proposal. 

Using instead of a maximal foliation near space-like infinity a double null 
foliation, Klaincrman and Nicolo proof a similar though technically simplified 
result (lElj). Revisiting their proof in the light of the results discussed above 
they give asymptotic conditions on the data near space-like infinity which allow 
them to verify the peeling behaviour ('|92j). 

More recently Lindblad and Rodnianski obtained a technically even more 
simplified global existence proof by using the Einstein vacuum equations in wave 
coordinates (IHZI)- They avoid the difficulties near space-like infinity by starting 
from Cauchy data which are exactly Schwarzschild near space-like infinity but 
they plan to give a proof based on weaker assumptions in a forthcoming article. 
It is not clear yet to what extent their method will allow them to gain precise 
control on the smoothness resp. peeling behaviour of the fields near null infinity. 

Any analysis of the evolution of the fields near space-like infinity will need 
to impose restrictions on the initial data to obtain control on the evolution and 
on the behaviour of the fields near null infinity. In the articles j67j and (68j has 
been developed a setting which allows one to analyse under suitable regularity 
conditions on the data at space-like infinity the evolution of the fields near the 
critical set where space-like infinity touches the null infinities (a notion made 
precise in |^) at all orders. As a result it is shown that even for (conformal) 
data of maximal smoothness the solution can develop at all orders logarithmic 
singularities at the critical set and consequently on null infinity (cf. j7()|). Fur- 
thermore, there is obtained for the first time a series of analytical expressions 
which relates a certain class of obstructions to the smoothness at null infinity 
to certain specific fall-off properties of the initial data. While the analysis is 
carried out in an algorithmic way, the complexity of certain expressions grows 
so quickly with increasing order that the possible existence of a further class of 
obstructions was left open. 

Recently the remaining case was studied by Valiente Kroon by using an 
algebraic computer program f |132| '). He did find further obstructions and re- 
markably, up to the order to which the calculation could be performed, there 
is now evidence that in the case of time reflexion symmetric data smoothness 
(resp. C*^) at null infinity requires the data to be asymptotically static (resp. 
asymptotically static up to an order p for a certain integer p = p{k) which still 
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needs to be determined). Note that this is much weaker than 'static in a neigh- 
bourhood of space-Hke infinity', in which case that smoothness at null infinity is 
easily shown. For more general data the situation is not so clear yet but 

asymptotic stationarity may play an important role. A more detailed discussion 
of the situation is given in [721 • 

Since the work referred to above relies on a conformal respresentation of the 
Einstein equations which only seems to be useful in 4 dimensions (cf. |69| ) . it 
may be worth mentioning here that other conformal field equations have been 
suggested recently by Anderson which work in all even space-time dimensions 
(0)- It can be expected that many of the results obtained in 4 dimensions 
can be generalized to all even dimensions Recent results by Hollands and 

Wald (IHll) and Rendall f lllTp suggest, however, that conformal equations with 
similar properties do not exist in odd space-time dimensions. 

Besides clarifying the asymptotic behaviour of gravitational fields one of the 
main motivations for the work in 68 was to provide a setting which would allow 
one to calculate numerically entire space-times on finite grids, including their 
asymptotic structure and radiation fields. Once the numerical evolution can 
been pushed past the critical set the solution will contain hyperboloidal slices 
and earlier analytical and numerical results can be applied. The numerical 
implementation for this program has not been given yet. As shown by the 
following example, there arc interesting question which would be difficult to 
study by other numerical approaches. 

Beig and O'Murchadha describe in J7] an interesting construction of asymp- 
totically flat initial data on 5 = R'^ with trapped surfaces and suggest that 
these surfaces are due to concentration of gravitational radiation (these data 
are not obtained by using initial hypersurfaces with non-trivial topology but by 
analysing sequences hn of conformal metrics with positive Yamabe number on 

for which the Yamabe number approaches zero as n — > oo). Since the data 
are time reflection symmetric, the singularity to be expected in the future must 
be considered as a reflection of the singularity in the past and cannot be inter- 
preted as being due to radiation. In 18 similar data are constructed without 
the time reflection symmetry, but this by itself docs not preclude the existence 
of a singularity in the past. 

Dafermos has recently shown the existence of maximal developments arising 
from asymptotically flat Cauchy data for the spherically symmetric Einstein- 
scalar field equations, which contain an event horizon in the future but for 
which all causal geodesies are complete in the past (^48^). One would consider 
a solution to the Einstein vacuum field equations with these global properties 
as presenting a black hole due to a collapse of gravitational radiations. Clearly, 
it is an interesting question whether such solutions do exist or whether they are 
excluded by the field equations. 

This question cannot be answered by analysing spherically symmetric situ- 
ations. Since the techniques used in 0H1 only apply to wave equations in two 
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space-time dimensions, the answer is not known. If one could perform numeri- 
cal calculations which cover the maximal globally hyperbolic solution space-time 
one might be able to show the existence of solutions as indicated above. A nat- 
ural problem to consider here is the characteristic initial value problem for the 
conformal vacuum field equations where data are prescribed on a cone repre- 
senting past null infinity (cf. |64j ) . In that case one would have perfect control 
on the past but there arise other difficulties (non-smoothness of the initial hy- 
persurface, how to prescribe the data on the cone to obtain an asymptotically 
flat solutions, development of caustics, the transition of the numerical evolution 
process through space-like infinity, etc.) which let the setting indicated above 
look more attractive. 

We end this article by pointing out a problem which did not receive much 
attention yet but which will become important as soon as certain technical ques- 
tions, which are still under investigation, will be understood. Depending on the 
way it is approached, it poses itself differently, but it is most severe if radia- 
tion fields are to be calculated numerically. There are essentially three different 
approaches to the numerical calculation of radiation fields: (i) the standard 
approach based on Cauchy data and the introduction of an artificial time-like 
boundary to make the computational grid finite, (ii) semi-global approaches 
based either on characteristic foliations and characteristic initial hypersurfaces 
extending to null infinity or on the conformal field equations and hyperboloidal 
hypersurfaces, (iii) global approaches, like the one indicated above, which as- 
pire to calculate entire space-times (possibly including their asymptotics). In 
all three cases inner boundaries may be considered to avoid the approach to 
singularities, but such boundaries will be ignored here. 

In the standard approach particular technical problems are introduced by 
the presence of the time- like boundary, which, being in general not distinguished 
geometrically, is somewhat unnatural. Nevertheless, it has been shown by Nagy 
and Friedrich that the vacuum field equations admit well posed initial-boundary 
value problems (which includes, of course, that all constraints be satisfied). 
There has been discussed the freedom to prescribe boundary conditions and 
data and the pricipal difficulties and specific features of the problem have been 
pointed out {^T^). 

It turns out that three real functions can be prescribed on the time-like 
boundary. In the setting considered in |74| these are the mean extrinsic curva- 
ture, which can be understood as controlling the evolution of the boundary in 
the solution space-time (it does not suffice to say 'the boundary is the hypersur- 
face {x = 0}' for some coordinate x), and the other two are components of the 
conformal Weyl tensor which may be interpreted as controlling the two radiative 
degrees of freedom, though in general a fully satisfactory physical interpretation 
of these data does not exist. 

As pointed out in scction l3.il a non-constant mean extrinsic curvature cre- 
ates problems in the analysis of the constraints on space-like hypersurfaces. 
Therefore is is worth mentioning that the choice of a non-constant mean ex- 
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trinsic curvature on the time-like boundary also creates certain difficulties (cf. 
\74]). Though these are quite different from the ones encountered on space- like 
hypersurfaces, it seems to indicate that something important about the mean 
extrinsic curvature is not understood yet. We shall not be concerned here with 
this question, however, because it does not prevent us from solving the initial- 
boundary value problem in all generality. 

While the representation of the field equations considered in has been 
used in numerical calculations before (cf. [HH), it is not the one used in the 
majority of general relativistic numerical codes and there is now a considerable 
amount of work being done to derive similar results based on other main evolu- 
tion equations (cf. |124| and the references given there). It can therefore 
be expected that the basic numerical problems arising from the initial-boundary 
value problem will soon be overcome. Moreover, numerical experiments will also 
show how the boundary data must be prescribed to ensure a regular long time 
evolution of the boundary. 

One will then have to provide a meaningful concept of outgoing radiation in 
an initial-boundary value problem. A well-defined rigorous definition is not in 
sight since in general there does not exist a distinguished outgoing null vector 
field transverse to the boundary (the spherically symmetric case is trivial and 
in no way representative). We shall not be concerned with this question here, 
pretending that some answer can be given. Then there will still remain the 
question: How should one dispose of the freedom to prescribe boundary data ? 
This problem has hardly been considered so far and physical intuition is not 
likely to give an answer. 

That I am not raising here a purely academic question is illustrated by some 
recent calculations. Allen et al. intend to study in the article |21 radiation 
tails for black hole evolutions by solving numerically initial boundary value 
problems for wave equations on a Schwarzschild background. The problem is 
readily reduced to a problem for a wave equation of the type dudyip = Ui(r) ip, 
where u and v denote the standard retarded and advanced Schwarzschild time 
coordinates. The authors prescribe certain initial data and impose the boundary 
condition 

d^i^^O on T = {r = ro}, (3.11) 

for some suitable value > 2 m of the standard Schwarzschild coordinate r. It 
turns out that radiation tails as predicted by Price's law ( 113^) cannot be iden- 
tified in the subsequent calculations. The authors conclude in their summary: 
We have shown that finite-radius boundary conditions prevent the formation of 
power-law tails. 

The results of these numerical calculations have been confirmed analytically. 
It has been shown by Dafermos and Rodnianski (|Sn|) that in the setting of |2] 
any tails vanish on the event horizon faster than p{v)^^ as w — > oo, where p is 
any polynomial of the advanced time coordinate v. 

One should not think, however, that these rigorous results confirm the con- 
clusion above. Gundlach et al. successfully verify Price's law numerically by 
solving characteristic initial value problems on the same background (IHIl). This 
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is in accordance with the analytical results of 0^1 > where the characteristic ini- 
tial value problem is analysed with data prescribed on a pair of intersecting 
null hypersurfaces Nout, Nin which extend to future null infinity and across the 
horizon respectively. 

It follows now from the general theory of initial-boundary value problems 
that the data given on Ni„ and the data induced on a time-like hypersurface T 
which extends from Noutf^Nm to future time- like infinity determine the solutions 
considered in |SJ and HHI uniquely in the future domain of dependence of 
T U Nin (the set of points p for which any inextendible past directed causal 
curve through p meets T U Nin). There is no reason why the solution induced 
on D+ by the numerical results of [ST] should not be reproducible by a numerical 
calculation based on the data on T and iVj„. For this purpose one would have 
to require the boundary condition 

dyilj = d on T, (3.12) 

where the function d on T would have to be read off from the solution given by 

m 

The calculations in "7 thus do not indicate that power-law tails cannot be 
calculated by solving initial-boundary value problems, they just confirm that 
there is a serious problem: In general, one does not know the 'correct ' boundary 
data on the right hand side of H3.12|l . 

It may be said that the radiation signals one is really interested in will not 
be as delicate as the radiation tails. But this does not tell us how strongly 
changes in the boundary data will affect wave forms and, in particular, what 
will happen in a long time calculation. By a wrong choice of boundary data, 
the system which was to be modelled may be affected so drastically that the 
approximate radiation field calculated in the end has little do with the system 
envisaged originally. To what extent this can be the case can probably only be 
explored by numerical experiments. 

The discussion above seems to indicate that the calculation of wave signals 
based on characteristic, hyperboloidal, or standard Cauchy problems will more 
robust than those based on initial boundary value problems. But what ever 
one does, there will always be an arbitrariness in the choice of data. The 
work initiated by Corvino, discussed in section ITTl clearly illustrates the large 
freedom to deform the data outside compact sets, or, in other words, how to 
choose the extension near space-like infinity of the data which characterize 'our 
system' on a compact set. 

This leaves one with the task of minimizing the import of accidental infor- 
mation. The notion of spurious radiation may have an intuitive meaning but it 
cannot be well defined. There are two analytical suggestions, however, which 
may supply useful criteria. Dain was able to associate with a given Cauchy data 
set a certain number, related to the global structure in a similar sense as the 
total mass but defined in quite a different way, which vanishes if and only if the 
data are stationary ([HS)- This number may thus be considered as a measure 
for the radiation content of the data. It remains to be seen, how changes in 
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the data are reflected in changes of this number and whether minimizing this 
quantity in suitable classes of data leads to applicable results. 

The second suggestion follows from the analysis of asymptotic smoothness 
properties. As discussed above, smoothness requirements at null infinity imply 
asymptotic conditions on the Cauchy data near space- like infinity like asymptotic 
staticity or asymptotic stationarity. These may be interpreted as suppressing 
'spurious radiation' near space-like infinity. 

In any case there will remain some freedom and one needs to assess the 
effect of changes of the data near space-like infinity on the wave forms. In 
the end one may have to look for features of radiation signals which are stable 
under the remaining admissible changes of the data. Identifying such features 
and understanding to what extent they characterize the type of the source is 
certainly an important task. It will require numerical as well as analytical input. 
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